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Abstract 

In this work, based on a realization of an inhomogeneous Poisson 
process whose intensity function depends on a real unknown parame- 
ter, we consider a simple hypothesis against a sequence of close (con- 
tiguous) alternatives. Under certain regularity conditions we obtain 
the power loss of the score test with respect to the Neyman-Pearson 
test. The power loss measures the performance of a second order ef- 
ficient test by the help of third order asymptotic properties of the 
problem under consideration. 
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1 Introduction 

Let X^ n > be a realization of a nonhomogeneous Poisson process observed on 
some increasing subsets A n , n = 1,2,... of d dimensional Euclidian space 
M. d with intensity function S (i?, x) , x 6 A n depending on one-dimensional 
parameter •§ G 0. Based on we want to test the hypotheses 

Hi : > #o, 
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where $0 is a given value in the parameter space 9 . Let us fix some a G (0,1' 
and define the class K,^ of tests at level 1 

/Ci") = {0 n : Etf„0 n (IW) = a) 



where E# denotes the mathematical expectation with respect to the proba- 
bility measure P^. With fixed n, generally speaking, there is no uniformly 
most powerful test in K,^ (see [8]). Therefore we turn to the asymptotic 
approach and introduce the class K' a of sequence of tests of asymptotic level 
1 — a, i.e., 

K! a = |{<M: lim E*A = a) . 

It is well known that if n — > 00 for any given value d of the alternative 
the power of any reasonable (consistent) test tends to 1 ( see [?], PQ). In 
order to compare the different tests we use the Pitman's approach (see 
[7]) where instead of a fixed alternative 1? we consider a sequence of so- 
called local alternatives (close or contiguous alternatives) which converges 
to $0 with a certain rate and hence it is difficult to distinguish between 
the null hypothesis and alternative. More precisely, let {<£>„} be a sequence 
of nonnegative numbers which converges to zero with such a rate that the 
likelihood ratio 

jp(n) 

Z n (u) = (IW) 

has a nondegenerate limit for any u with i?q + ip n ii G 6. By the Neyman- 
Pearson lemma the most powerful test for Ho : i? = $0 against the local 
alternative 7i u : $ = $0 + y? n -u with w > 0, is given by 



0n (I W ) 



1, if A n (u) > b n (u) 
0, if A n (u) < b n {u) 



where A n (w) = h\Z n {u) and the constant b n {u) together with the contribu- 
tion of the randomized part provide the size E$ (fi n (X^) = a. The power 
of (fi n as a function of u is called the envelope power function . For any fixed 
n it is the supremum of the power at the local alternative d u = $0 + <£nU 
over all the tests at level 1 — a, i.e. 

E#J n (iM) = sup E^<fi n (iW) . 

Notice that (fi n is not a test for the main hypotheses Ho and Hi because it 
depends on the parameter u. 
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Let us introduce the score statistic 



where the normalizing factor ip n = J„( , (9o) _1//2 is the inverse square root of 
the Fisher information 

T (A) - f g ^' X ) 2 dr 

Wo) " J An S(* ,z) dX 

at the point iV Here ir^ (dx) = (dx) — S (i?q, x) dx is the centered 
Poisson process (for the definition of a stochastic integral w.r.t. a Poisson 
process see the next section) and S x) denotes the derivative of S x) 
with respect to f). Based on A n ($ ) we introduce the score test 



1, if A n (t? ) > z 
0, if A n (0 O ) < z 



where z a is 1 — a quantile of standard Gaussian law, i.e., P {( > z a } = a 

and C ~ A/"(0, 1). It is well known that if the family |p^ n) , e ej of 

distributions is locally asymptotically normal (LAN) at the point $0; then the 
test 4> n e JC' a is locally asymptotically uniformly most powerful (LAUMP) 
(or first order efficient), i.e., for any K > 



sup 

0<u<K 



E*J n (XW)-E<>J n (xW) =o(l) 



as n — > oo (see |12j). Moreover the power function of <ft n at admits the 
representation 

E*J n (I W ) = P {C > 2a - «} + O (1) (1) 

for any u > 0, where £ ~ A/"(0, 1). For n large, hence the power of (j) n 
approximates the envelope power function up to order o(l). A refinement of 

(0Q) is given in PJ. For a family jpjf \ i? E O j of distributions related to a 

Poisson process with intensity functions {S •) , i? G 6} the conditions of 
LAN for multidimensional parameter d are obtained by Yu. A. Kutoyants, 
[B]. The first order efficiency of 0„ follows from the LAN representation 
which implies in turn the asymptotic normalities of A n (i? ) and A n (u) under 
both Tfo and 7f u . Therefore the refinement of the central limit theorem, by 
taking into account one term after the Gaussian term, improves the situation. 
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This can be done by the help of Edgeworth type expansion of the distribution 
function of the stochastic integral A n ($ ) and A n (u) under TCo and TC U . Under 
certain regularity conditions related to second order asymptotic properties 
of the family jp£°, G ej, we can construct a second order efficient test, 
i.e., a test 0* such that for any K > 

= o(4), (2) 

for some sequence e n — > (see [1]). Furthermore the probability of the 
first type error of 0* is given by E# o 0* (X (n) ) = a + O(e^). From ([2D it 
follows that the power function of 0* (a second order efficient test, generally) 
approximates the envelope power function up to order 0(e^) and hence it 
works as good as 0„ up to this order. Second order efficiency of 0* is related 
to the fact that the two first terms in the Edgeworth expansions of the 
distributions functions of A^^o) and A n (it) under the local alternative are 
equal up to the order 0(el). Hence to measure the performance of a second 
order efficient test, and especially 0* , it is natural to consider the power loss 
of 0* with respect to the most powerful test n , which is defined by 

r(«) = lim e- 2 (e*A (*<»>) - E,. u 0; (X^)) , (3) 

for u > 0. This requires to take into account higher order terms in the 
Edgeworth expansions of the distribution functions. See pQ, chapter 3, for 
a general theorem and the power loss results for the tests based on L— , R— 
and U— statistics in the i.i.d. case. The main object of this work is to 
obtain the power loss of the score test 0* based on a realization X^ n ' of a 
nonhomogeneous Poisson process with intensity function S (#, x) , x G A n 
and to give the explicit representation of r(u). 



sup 



E, u n (JW) - E^0; (IW) 



2 Preliminaries 

Let us remind several facts from a Poisson process. A Poisson process 
is a random point measure which on the set B C A n can be written as 

where {x,} are the events (random points) of the Poisson process and x^ B y is 
the indicator function of the event D. The Poisson process with (parametric) 
intensity function S($,x), x G A n (with respect to Lebesgue measure) is 
entirely defined by the following two conditions: 



4 



• for any collection of disjoint sets Bi, . . . , B m C A n the random variables 
X( n )(Bx), . . .,l( n )(B ffl ) are independent, 

• the random variable X^(M) for any B C A„ has Poisson distribution 
with parameter A^(B) = J M S(fi, x) dx. 

This and other definitions of the spatial Poisson processes as well as their 
properties and examples can be found in many books devoted to point pro- 
cesses (see, e.g., Daley and Vere- Jones (2 . Krickeberg [SJ, Reiss [ID] . Ripley 
[11] . Snyder and Miller [IS])- The spatial Poisson processes are widely used 
in many fields. In particular this is a good mathematical model for the prob- 
lems of image restoration when the optical signal is weak and statistics of 
photons is well described by an inhomogeneous Poisson process [13] . 

By the definition, is a random element of the set A4q , containing all 
integer valued measures defined on the set A„. Let *S(J\Aq ) be the smallest 
a— field with respect to which all the mappings: 

n B : - {0, 1, 2, ... , oo}, n B (lW) = lW(B), B G A n , 

are measurable. Here A n is the a— field of Borel subsets of A n . Let Pj^ de- 
note the probability law induced by the random element (realization) of 
a Poisson process with intensity function S(i9, x), x G A n on the measurable 
space (M^\^8(M^)). We remind that if the intensity measures A£ } and 
A^ are equivalent then the corresponding probability measures and 

(n) 

P$ 2 are equivalent and the likelihood ratio is given by 

( X (n)\ _ 

= exp(/ ln^^-X^(dx)- [ [S($ 2 ,x)-S(A,x)]dx 
lJA n &(yx,x) J An 

For a proof see [6] page 28. Here the stochastic integral 

f f(x)X^(dx)= f( x i) 

where {x{\ are the events (random points) of the Poisson process. 

Since the main tool used in this work is based on the Edgeworth expan- 
sion, here we present the conditions under which the distribution function 

F n {y) = vf {/„(/) <y}, 
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of the stochastic integral 



/»(/)= / fn(x) 7rW(dx), 



admits an Edgeworth type expansion for two terms after the Gaussian term, 
where ir^ (dx) = 1^ (dx) — S($, x) dx is the centered Poisson process. We 
suppose that 



/ f n (x) 2 S($,x)dx = l. 

J An 



The expansion is obtained under the following two conditions: 

B\. There exists a sequence of real numbers e n — > 0, as n — > oo and con- 
stants C r > 0, r = 3, 4, 5, such that 



/ 

J A, 



|/ n (x)| r S(tf,x) dx<C r e r n 



r-2 



B 2 . There exist constants 7 > 5/2 and c > satisfying the inequality 
ff c + ffc 2 , + eg - \ < such that 

inf / sin 2 (t/ n (x)) 5(7?, x) da; > Tine' 1 

2 2 

for all large n. 
Let us introduce the cumulants: 



lr,n = / fn{x) r S(0,x) dx, 



r = 3,4 



and the Hermit polynomials: 

H 2 (y) = y 2 -l, H 3 (y) = y 3 - 3y, H 5 (y) = y 5 - 10y 3 + 15y. 
Theorem 1. Let the conditions Bi, B2 be fulfilled, then uniformly in y G 



2 

^(v) = my) - |f " if Fs(y) _ I2 HM n{y) + 0(4) ' 

for all n large. Here Af(y) and n(y) denote the distribution and density 
functions of the standard Gaussian law, respectively. 
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Note that 7 r , n = 0(e r n ~ 2 ), r = 3,4 by B x . 

Proof. See [3], Page 36. The proof is a special case of a general theorem 
given by Kutoyants, where the expansion is obtained by the powers of e n up 
to order e*, k = 1,2,... (see ([6]), page 131). 

Corollary 1. Let < a < 1 be given and the conditions B\ — B2 be fulfilled. 
Then the equation F n (y) = 1 — a + 0(e n ) has a solution y = c n ^ a , 

2 

C n ,a = Z a + ^(^a) + ""777"" ^(^a) + "St -^sOO- 

o 24 72 

See [3], page 40 for proof. 



3 Second Order Efficiency 

Let j3 n (u, 4>n) denote the power of a test <p n at the local alternative $ u = 
$0 + VnU, i.e., 

Under certain regularity conditions, slightly weaker than T>\ — D 3 (see the 
next section and |4J), the score test 



based on 



with 



b * ( X (n)\ = f 1 ' ifAn(^o)>c n 

n[ ' \0, if A n (i? ) < Cn, 



73,n /, 2\ 3 /* g (^0^) 3 , 
Cn = Z a — (1 - Z a ) , 7 3jn = (f n I dx 
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is second order efficient, i.e., it satisfies 



sup 

0<u<K 



for any K > 0. For proof see pE], Theorem 6. Indeed to establish the 
second order efficiency we obtain the following representation of distribution 
functions of A„ ($0) and A n (u) under the local alternative d u : 

Pn [U, <P n ) = M ( '- (1 - {U - Z a ) ) n(u - Z a ) + 0{S n ) 



Vn 



Pn(u, (pn) = AM pr (1 - [U - Z a ) )n(u - Z a ) + 0{E n ), 

V aJu) I 6 



where 



m n (u) = E, u A n (0 O ) = <p n [ ( S W»> X ^~ S *)) dx, 

i£ = E, u (A n (tf ) - m n {u) f = &( (0 U , x) dx, 

Ja„ S (tf , x) 

» n (u) = E# u A n (u) = [ fin (0„, x) - S (0 tt , x) + 5 (0 O , x)) dx, 

Ja„ V ^ (yo, x) j 

o* n (u) = E*. (A n (u) - ^(u)) 2 = / fin 2 5 (0 tt , x) dx, 

73,n(«) = -3 / (tf u ,x) dx, 



Hence it suffices to show that 



sup 1 7 3 ,n (u)--f' 3 (u)\=0 (e 2 n ) , 

0<u<K 



sup 

0<u<K 



m n (u)-c n n n {u)-b n (u) 



Vn <T n (u) 



= 0{el). 



The case u = corresponds to the size of <p n which is equal to a and the size 
of 0* given by a + 0(4). The latter follows from the Edgeworth expansion 
of distribution function of A n ($ ) under the null hypothesis. 

Representation of the power . Now we obtain the explicit represen- 
tation of the power of 0* (up to order 0(4))- The Taylor expansion 

2 2 3 3 

S(# u ,x) = S(# ,x) + <p n uS(# ,x) + ^-S^ ,x) + ^Ls^($ n ,x), 
for some intermediate point i?o < $n = $n(u, x) < implies that 

ro „ (M) _ c „ = „ _ Za + £ s( ^g^> d, + ^(i - 4) + (4). 

Similarly we obtain 

^ = 1 + jf dx + 0(4) = 1 + M73 ,„ + 0(4). 
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Hence 



m n (u) -c n 1 - zl - 3u (u - z a ) 
= u - z a + 7 3jn + 

Vn 6 

2 J An S(# ,x) dX + U 



n J 



On the other hand by the help of Taylor expansion we obtain 

73,n(w) =l3,n + 0(e 2 n ). 

Therefore we obtain the following representation 

(3 n (u, <p* n ) =M{u-z a ) + Q n (u) n{u - z a ) + 0(4), (4) 
for any u > 0, where the polynomial (in u) 

u(z a -2u) ip\u 2 f S(4i)5^o,i) j 

Qniu) = 7 3 ,„ + — — / -T- r dx, 

is of order 0(e n ). The equation (jlj), refines the first order representation (JT]). 
Notice also that the the second order efficiency of 0* implies that we have 
the same representation for the power of n . 



4 Power Loss 

The power function of a second order efficient test agrees with that of the 
most powerful test up to order O(e^). Hence it is natural to consider the 
power loss of </>*, which is defined for any u > by 

r(u) = lim e~ 2 (f3 n (u, 0„) - /3 n {u, fa)) . 

n— >oo \ / 

Below S^'^jX) denotes the jth derivative of S^, x) with respect to We 
write x) and S($, x) for the first and second derivatives, respectively. 
We consider the following conditions: 

T>\. The intensity function S{$,x) is four times differentiable with respect 
to i? in a right neighborhood of ^o- 

T>2- The conditions B\ and B2 are satisfied for the stochastic integrals 
A n ($ ) and A n (u) under Ti and Ti u with some sequence e n — > 0, 
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V 3 . There exists some functions fj(x),j = 0,...,4, x G A n not depending 
on such that S(0,x) > f (x), \S^(0,x)\ < fc(x), j = 0, 4 for all 
x G A„ and all in a right neighborhood of O . We suppose also that 

, /■ l/i(s) |fc 



Example 1. Let be a realization of a Poisson process on the set 

A n = [0,n] with positive intensity function S(0, x) = 0S(x) + A (amplitude 
parameter) or S(0, x) = S(0 + x) + A (phase parameter), where S(-) is a two 
times differentiable periodic function and A > (dark current) is a known 
constant. In both cases the condition T> 3 is satisfied with ip n ~ C rT 1 ! 2 
for some C > and e n = n^ 1 ' 2 . For the frequency modulation model 
S(0,x) = S(0x) + A we have ip n ~ C n~ 3/2 and e n = n~ 1/2 . 
Let J n denote the quantity: 

's (0 o ,xf- S (0 o ,x) S (0 o ,x)) 2 

K = J —dx- 

S(0 o ,xf 

S (0 O , x) (s (0 O , xf - S (0 O , x) S (do, x)) 
S(0 O , x f 

Note that J n = O(e^), by T> 3 . We have the following theorem: 

Theorem 2. Let the conditions T>\ — T> 3 be fulfilled. Then the power loss of 
0* with respect to the most powerful test <p n is equal to 

r(u) = J.), 

for any u > 0. 

Proof. By T> 2 we can write the following third order expansions: 
f3 n («, <Pl) = M(a n ) + H 2 (a n ) n(a n ) - H 3 (a n ) n(a n )- 
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H 5 (a n )n(a n ) +0{e 6 n ) 



P n (u, 4> n ) = Af(A n ) + H 2 (A n ) n(A n ) - H 3 (A n ) n(A n )- 



7 3 ,nH 

72 



2 



^ 5 (A n )n(A,) + 0(eJ[), 
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where 



_ m n {u) - c n _n n (u)-b n (u) 



= ^Mr^ S (^,x)dx, r = 3,4 

Vn JA n D {V , X) 

Using the Taylor expansions of S ($„, x) and In fp^fy, we get 

74,n(«) = 74,n + 0(4), 73,n(«) 2 = Ta.n + °(4) 
74,n(«) = 74,n + 0(4), 7 3,» 2 = 73,n + 0(4)- 

Since A n — a n = 0(4) (which follows from the second order efficiency of </>*), 
then 

M (AO = M K) + (A, - On) n(a n ) + 0(e 4 n ). 

Therefore 

Pn (u, <t>* n ) = N(a n ) + ^jf H 2 (a n ) n(a n ) - ^ H 3 (a n ) n(a n )- 



2 

H 5 (a n )n(a n ) + 0(e 3 n ), 



72 

73 n ( M ) 

Pn(u,4>n) = N {a n ) + (A n - a n ) n(a n ) + -^| — H 2 (a n )n(a n )- 

2 

- H 3 (a n ) n(a n ) - H 5 (a n ) n(a n ) + 0(4). 

Letting A = A (it) = u — z a and taking into account the fact that a n — A = 
0(e n ), we can write 

/5n(«, 4>n) - Pn («, O = (A» ~ a n) n(A) + 

+ g,(AWA) (7 > n(M) _ 73j , M)+0(4) . 

Therefore we have to consider the terms A n — a n and 73 n (u) — 73 in (u). Since 
the investigation is at accuracy level 0(4), we modify the threshold c n as 
follows: 

2 

c n = z a + H 2 {z a ) + H 3 (z a ) + H 5 (z a ). 
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With this constant the probability of error of the first kind of 0* is equal to 
a + 0(4)- Similarly the explicit form of b n (u), the threshold of n , can be 
written as: 

b n {u) =^ n + °n (z a + ^ H 2 (z a ) + ^ H 3 (z a ) + ^ H 5 (z a )\ + 0(4), 



where 



H n = E# A n (u) 



Ut X ) ■. 

J A n I S (&o,x) S (t? 0> a?) 



S (i? ) x ) da;, 



a 2 n = E# (A n (u) - n r , 



In 



Ul x , 



S (i? , dx, 



7r 



— / (\n S ^ u,x) ) S($ ,x)dx, r = 3,4. 



Below for the sake of simplicity we use the following notations: 

S($ ,x) r >S($ ,x) r2 



I(r ,r 1 ,r 2 ) 



S ($ ,x 



r 



-dx 



for nonnegative integers ro,ri,r%. The long but straightforward calculations 
show that (see Appendix) 



u 



2u 3 - 2u 2 z a - u(l - z 2 
9u 3 - I2u 2 z a - 6u(l - zl) 4 



4, -0„ = ^-^7(1,0,2)- 

+ 



a) A 



^/(2,2,1)+ 



24 



■¥>*/(3,4,0) + 



(5) 



9u*-Qu 2 z a + 2u(l-z 2 a ) 6 2 

+ 24 ^™ ^ ' ' > + 

6u 3 - 6u 2 z a - 3u(l - 4) 

-^^/(i,i,i) 2 + o(4)- 



+ 



12 



^ 1(2,3, 0)7(1,1,1)- 



Also for the term 73 (u) — 7 3i „ (w) we find that 

1?/ i7/ 1?/ 

73,n («) - 73,n («) = y ^ /(2, 2, 1) + - ^ 7(2, 3, 0) 2 - - <p' n 7(3, 4, 0)- 



37; 

- y ^7(2, 3, 0)7(1, 1,1) + 0(4). 



(6) 
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See Appendix for some details. Combining these results yield 

n~\A) 4> n ) - (3 n («, 0;)) = ^ ^ /(1, 0, 2) - ^ <p\ 7(2, 2, 1)+ 

+ y ^ A3, 4, 0) - ± ipi 7(2, 3, 0) 2 - £ ^ /(1, 1, 1) 2 + 
+ T ^ 7(2, 3, 0) 7(1, 1, 1) + 0(4) = - <4 [ 7(1, 0, 2) - 2 7(2, 2, 1) + 
+ 7(3, 4, 0)] - £ ^ [7(2, 3, 0) - 7(1, 1, 1)] 2 + 0(4). 



Hence for any w > 



(3„(u,(f> n ) - (3 n (w,0*) 



u 3 n(A) 
8 



7„ + 0(4) 



where 



J n =ti\ 7(1, 0, 2) - 2 7(2, 2, 1) + 7(3, 4, 0)] - <p 6 n [7(2, 3, 0) - 7(1, 1, l)] 2 . 
This can be simplified as 



T 4 



■dx- 



S(# , xf -S(# ,x)S(#o, x)J 
S(# ,x) (s($ ,x) 2 -S($ ,x)S($ , 



x , 



Note that by the Cauchy-Schwartz inequality; 

, S (#„, x) (S (#„, - ^ (#0, x) 5 (#0, 



dx 



dx < 



S(A,x) 



I 



-dx 



S(tf ,x) 2 -S($o,x)S($ ,x) 



/ 



(s($ ,x) 2 -S($o,x)S($ , 



S(0 o ,x) 



s(# ,xy 



-dx. 



-dx 



J 



Hence as one expect J n > 0, because the test </> n has the maximum of power 
at the local alternative Furthermore V3 implies that the quantity J n is 



13 



of order O(e^). Therefore 



(3 n (u,<j) n ) - Pn («,</>* 



u 3 n(A) 



e n J n + 0(e n ) 



which completes the proof of the theorem. 



Example 2. (Amplitude Parameter.) Suppose that we observe a 
realization of a Poisson process on the set A n = [0, n], n = 1,2, ••• 
with the intensity function 



S(ti,x) = $S{x) + \, 



> 



where A is a known positive constant (dark- current) and S(x) is a known, 
nonconstant, differentiate le (with respect to x) and periodic function with 
period r > 0. We have two hypotheses TLq : ■& = i? against 7Yi : i? > i?o> 
where $ > 0. The intensity function S($,x) is supposed to be positive in a 
right neighborhood of $o an d all x. We obtain the power loss of the test 



0* (X( n >) = X{A n (ti )>c n } 



with 



where 



2 

Z a + "~7T~ #2(^0) + -777 1 #3(0 + H 5 {z a ) 
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24 



72 



73, n 

74, n 

A 2 



S(x) 



1 T S(x) 4 



dx 



rA 4 Wo (tf S(x) + A) 3 
1 f T Six) 2 

ax. 



dx 



rj $ S(x) + \ 



The conditions T>\ — P 3 are satisfied (see [I], pp. 202-204) with e n = n l / 2 
and if n = rC 1 ! 2 A~ l {\ + 0(n -1 )) — > 0. One can obtain 



/ 



lim (e n 2 J n ) = r 



s(xY 



(t> S(i)+A) i 



dx 



JO {tf S(x)+\) A U - 



S(x) 2 , 



S(x) 2 



dx 



B. 



Hence the power loss is equal to 



u 3 n(u - z a ) D 
r(w) = fi. 
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Example 3. (Frequency parameter) In this example we consider a 
strongly nonhomogeneous case with nonclassical rate n~ 3 ^ 2 (instead of n~ l l 2 
in the i.i.d. case). Suppose that we observe a realization of a Poisson 
process on the set A n = [0, n], n — 1, 2, • • • with periodic intensity function 

S($,x) =e sin ^ x \ tf>0. 

We have the simple hypothesis TCq : d = $0 against Hi : 1? > ?V We consider 
the test 0; (JfW) 

= X{A„(^ )>cn} based on the statistic 
A n (0„) =<p n x cos($ x) (da;) - e sin(,?ox) da;) 
and the threshold 

2 

Cn = Z a + ^ #2^) + -777 1 ^(^a) + "If ^5(^0), 

o 24 rz 

where letting r = |^ we have 

cos 3 ($ x)e sin ^° x) da; 





c 3 


73,n = 






c 4 /• 


74,n = 




cr 2 = 


1 r 




3^7o 



cos 



! (t? x) e s ' m( ^ x) dx 



cos 



(^ a;)e sin(,9ox) da;. 



The conditions T>i — V 3 are satisfied (see [I], pp. 205-207) with e n = n 
and if n ~ C n~ 3 / 2 . 

The power loss is equal to 

9ru 3 n(u-z a ) C sin 2 (tf x) e sin ^ x ^ dx 
r{u) = — 



-1/2 



40 (£cos 2 (^ x)e si ^x) dx ) 2 ' 

for any u > 0. 

Representation of the power. We consider the explicit representation 
of the power of 0* up to 0(e 3 ). As we saw 

P n (u, <f)* n ) = N(a n ) + 2|p # 2 (a„) n(a n ) - -^p ffs(on) n(an) - 



2 

# 5 (a n )n(a n ) + 0(e 3 n ). 
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Since a n — A = 0(e n ), it can be written as 

P n («, 0;) = N{a n ) + ^ H 2 (a n ) n(a n ) - 77 3 (A) n(A)- 



7 77 5 (A)n(A) + 0(4). 
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Remind that A = u — z a . Indeed we have (see Appendix) 

an-A = y ^/(l,l,l) + y^/(l,l,0,l)+ 

1 — zi — 3u A o T .„ „ 
+ ^ <4 7(2, 3,0)+ 

9^ 2 A-2 M (1-^) ^ u ** in t T e> 9 n 
+ 24 V^n^(2,3,0) — </? n 7(2,2,l)- 

- T ^/(2,3,0)/(1,1,1) + O(e»). 

Taylor expansion yields 

N(a n ) = Af(A) + (an - A)n(A) - ^"^^ + 0(4) 

and 

7 3in 77 2 (a„)n(a n ) = 7 3jn 77 2 (A) n(A)+ 

+ A (2 - 77 2 (A)) n(A) 73 ,n K - A) + 0(4). 

Note that 

K " A) 2 = ^ /(1, 1,1)+ 1 -^~ 3mA v l 1(2, 3, 0)) 2 + 0(4) 
and 

73,n («n - A) = U - <p* n 7(2, 3, 0) 7(1, 1, 1) + <p* 7(2, 3, 0) 2 +O(4) 

Therefore we arrive at the following form 

(3 n (u, 0;) = M(A) + n(A) (n(n) + r 2 (n)) + 0(4), 
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where the terms r$(n) = 0(ejj, i = 1,2, 

n(n) = y ^ 1(1, 1, 1) + 7 3,n 

, . w 3 a . , 9w 2 A - 2w(l - z 2 ) 2 
r 2 (n) = - ^ 1(1, 1, 0, 1) + ^ 7 | n _ 

- — ipi 1(2, 2,1)-— ^ 7 3,n 1(1, 1, 1)- 

" y( y^/(l,l,l) + ^ 73, J + 

A(2-/7 2 (A)) /n 2 3 , l-* 2 -3^A _ 2 

+ g I y <^n 73,n i (1, 1, 1) + g 7 3 ,n 

2 

-^Jf 3 (A)n(A)-^# 5 (A)n(A). 
For the notations 7 above, see (J7J) in the Appendix. Especially 

\3 



^1(2,3,0) = ^ / §^4dx = 7 3,n. 



'a„ S (&o,x) 

Note also that ri(n) = Q n (u) given in the equation 



5 Appendix 

In order to obtain (jSJ) and (jSJ) we use the expansion of 5(i? u , x) about $0 an d 
the following Taylor expansions which precise the terms up to order 0(e 3 ). 
Letting 

_ S{ti u ,x) - S($ ,x) 

we have 

ln(l + t) = t - \t 2 + ^ 3 - -/ + O (t 5 ) 
2 3 4 v y 

ln 2 (l + t)=t 2 -t 3 + lV + 0(t 5 ) 

ln 3 (l + t)=t 3 -i 4 + 0(t 5 ) 
ln 4 (l + t) =t 4 + 0(t 5 ) . 
Note that the conditions T> 3 imply that for the terms of order (t 5 ) we have 

/ o (t 5 ) dx = o(4). 
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We use also the binomial expansions 

(1 + *)* = 1 + \s ~ \s 2 + O (s 3 ) 
(l + S )- = l-^ + ^ 2 + 0( s 3 ) 

for different values of s = s n (u) = 0(e n ) related to a\(u),a\ and rj^. We 
introduce the following notations: 



S(# ,x) ri S(# ,x) r2 (& ,x) 
S(0 o ,x) 



I( r o,n,r 2 ,r 3 ) = I p/Q ' da; (7) 



for nonnegative integers r ,ri,r 2 and r 3 G {0,1}. When r 3 = we write 
I(ro, n, r 2 ) instead of 7(r , n, r 2 , 0). 
Remind that 

H n (u) - b n {u) m n (u) - c n 



a n (u) r] n 

These give: 

2 4 3 4 
11 II II II 

M = - + - <p' n 7(1, 0,2) + y <pl 7(1, i,i) + -^ /(i, i, o, 1)- 

3 4 4 
II II II 

- - V l 1(2, 3,0) - — y, 4 /(2, 2,1) + - ^ /(3, 4, 0) + 0(4) 
= u" 1 - | ^ /(l, 0, 2) - I y, 3 /(l, 1, 1) - £ ^ /(1, 1, 0, 1) + 
+ 2 i ^/(3,4,0) + T ^/(l,l,l) 2 + O(4) 

W = — + ^ /(l, 0, 2) - ^ 7(1, 1, 1)- 

■u 4 — -u 3 2;„ , w , 2-u 3 — w(l — 2 2 ) — 3w 2 z Q , w 
g— ^ /(1, 1, 0, 1) + tfl A2, 3, 0) + 

2^- t i»(l-^)-3u 3 ^ 4 
+ ^ Vn J ( 2 > 2 > !)- 

6m 4 - 6m 2 (1 - z 2 a ) - 11m 3 z a 4 
24 V?„ J (3, 4, 0)- 

4M 2 (l-z 2 )+3 M 3 z a 6 2 

¥ 3 n-'( 2 >3,0) + 



24 

+ ^"f^ 3 " 3 ^ <P* n 1(2, 3, 0) 7(1, 1,1)-^ V l 7(1, 1, 1) 2 + 0(. 3 ). 
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Hence we arrive at 

7/3 -7#2 y3 

A n = u-z a + -^ n 1(1, 0, 2) + - ipl 7(1, 1, 1) + - ^„ 1(1, 1, 0, 1)- 
3n 2 -(l-a-_3^^ /(2)3)0) _ 



6 

<4 7(3, 4,0)+ 



3u 3 - 6m(1 - z 2 a ) - \2u 2 z a , 



24 

+ Mi-.~W % gj(23o)2+ 

+ 3tt '- 3u(1 ;;' ) - 6u ^ rf/(2,3,o)/(i,i,i)- 

-^rf/(i,i,i) 2 + o(4). 

For the term a n we have 

7/2 ^,3 

m» = « + — <p' n 1(1, 1, 1) + - 7(1, 1,0,1) + 0(4), 

C 1 = 1 - I ^ A2, 3, 0) - ^ ^ 7(2, 2, 1) + ^ tpl 7(2, 3, 0) 2 + 0(4). 

Therefore we can write 

a n = u - z a + — <p 3 n 1(1, lA) + -r <P$ 1(1, 1, 0, 1)+ 
z o 

+ ^-g^ 1(2, 3, 0) + 

9n 2 (w - - 2m(1 - zp 6 2 n 2 (?/-z a ) 4 
H ^ ^^(2,3,0) ^7(2,2,1)- 



u 3 



24 



- T ^ 7(2,3, 0)7(1, 1,1) + 0(4 

Hence one can obtain A n —a n as in (jHJ). Using the necessary Taylor expansions 
we get 

3m 

73 ,n(«) = 7(2, 3, 0) + « ^ 7(3, 4,0) - y ^, 7(2, 3, 0) 2 + 0(4), 
j' 3jn (u) = ipl 7(2, 3,0)-- 7(3, 4,0) + y 7(2, 2, 1)- 
-y ^7(2, 3, 0)7(1, 1,1) + 0(4), 
which prove ([6]). 
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